A CMV-Based Eigensolver for Companion Matrices

نویسندگان

  • Roberto Bevilacqua
  • Gianna M. Del Corso
  • Luca Gemignani
چکیده

In this paper we present a novel matrix method for polynomial rootfinding. We approximate the roots by computing the eigenvalues of a permuted version of the companion matrix associated with the polynomial. This form, referred to as a lower staircase form of the companion matrix in the literature, has a block upper Hessenberg shape with possibly nonsquare subdiagonal blocks. It is shown that this form is well suited to the application of the QR eigenvalue algorithm. In particular, each matrix generated under this iteration is block upper Hessenberg and, moreover, all its submatrices located in a specified upper triangular portion are of rank two at most, with entries represented by means of four given vectors. By exploiting these properties we design a fast and computationally simple structured QR iteration which computes the eigenvalues of a companion matrix of size n in lower staircase form using O(n2) flops and O(n) memory storage. So far, this iteration is theoretically faster than the fastest variant of the QR iteration for companion matrices in customary Hessenberg form. Numerical experiments show efficiency and accuracy of the proposed approach.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Lightweight 4x4 MDS Matrices for Hardware-Oriented Cryptographic Primitives

Linear diffusion layer is an important part of lightweight block ciphers and hash functions. This paper presents an efficient class of lightweight 4x4 MDS matrices such that the implementation cost of them and their corresponding inverses are equal. The main target of the paper is hardware oriented cryptographic primitives and the implementation cost is measured in terms of the required number ...

متن کامل

Compression of unitary rank-structured matrices to CMV-like shape with an application to polynomial rootfinding

This paper is concerned with the reduction of a unitary matrix U to CMV-like shape. A Lanczos–type algorithm is presented which carries out the reduction by computing the block tridiagonal form of the Hermitian part of U , i.e., of the matrix U +UH . By elaborating on the Lanczos approach we also propose an alternative algorithm using elementary matrices which is numerically stable. If U is ran...

متن کامل

Toward High Performance Divide and Conquer Eigensolver for Dense Symmetric Matrices

This paper presents a high performance eigensolver for dense symmetric matrices on multicore architectures. Based on the well-known divide and conquer (D&C) methodology introduced by Cuppen, this algorithm computes all the eigenvalues of the symmetric matrix. The general D&C can be expressed in three stages: (1) Partitioning into subproblems, (2) Computing the solution of the subproblems and (3...

متن کامل

Darboux transformations for CMV matrices

We develop a theory of Darboux transformations for CMV matrices, canonical representations of the unitary operators. In perfect analogy with their self-adjoint version – the Darboux transformations of Jacobi matrices – they are equivalent to Laurent polynomial modifications of the underlying measures. We address other questions which emphasize the similarities between Darboux transformations fo...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 36  شماره 

صفحات  -

تاریخ انتشار 2015